The garden of forking paths…

Currently stuck proliferating options.

Covariate options

## Loading required package: lattice
## Loading required package: latticeExtra
## Loading required package: RColorBrewer
## 
## Attaching package: 'latticeExtra'
## The following object is masked from 'package:ggplot2':
## 
##     layer

Resulting data

## ── Attaching packages ─────────────────────────────────────────────────────────────────────────────────── tidyverse 1.2.1 ──
## ✔ tibble  2.1.3     ✔ purrr   0.3.2
## ✔ tidyr   1.0.0     ✔ stringr 1.4.0
## ✔ readr   1.3.1     ✔ forcats 0.4.0
## ── Conflicts ────────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ purrr::accumulate()     masks foreach::accumulate()
## ✖ tibble::as_data_frame() masks dplyr::as_data_frame(), igraph::as_data_frame()
## ✖ purrr::compose()        masks igraph::compose()
## ✖ tidyr::crossing()       masks igraph::crossing()
## ✖ tidyr::expand()         masks Matrix::expand()
## ✖ tidyr::extract()        masks raster::extract()
## ✖ dplyr::filter()         masks stats::filter()
## ✖ dplyr::groups()         masks igraph::groups()
## ✖ dplyr::lag()            masks stats::lag()
## ✖ latticeExtra::layer()   masks ggplot2::layer()
## ✖ tidyr::pack()           masks Matrix::pack()
## ✖ .GlobalEnv::select()    masks dplyr::select(), raster::select()
## ✖ purrr::simplify()       masks igraph::simplify()
## ✖ tidyr::unpack()         masks Matrix::unpack()
## ✖ purrr::when()           masks foreach::when()
## 
## Attaching package: 'lubridate'
## The following object is masked from 'package:igraph':
## 
##     %--%
## The following object is masked from 'package:base':
## 
##     date
## OGR data source with driver: ESRI Shapefile 
## Source: "/Users/mrajeev/Documents/Projects/MadaAccess/data/processed/shapefiles/mada_communes.shp", layer: "mada_communes"
## with 1579 features
## It has 26 fields
## OGR data source with driver: ESRI Shapefile 
## Source: "/Users/mrajeev/Documents/Projects/MadaAccess/data/processed/shapefiles/mada_districts.shp", layer: "mada_districts"
## with 114 features
## It has 23 fields
## Joining, by = c("date_reported", "id_ctar")
## Joining, by = "id_ctar"
## Joining, by = c("year", "id_ctar")
## Warning: Column `distcode` joining factors with different levels, coercing
## to character vector
## Warning: Column `group_name`/`distcode` joining factor and character
## vector, coercing into character vector

## Warning: Column `group_name`/`distcode` joining factor and character
## vector, coercing into character vector

## Warning: Column `group_name`/`distcode` joining factor and character
## vector, coercing into character vector
## Warning: Column `commcode` joining factors with different levels, coercing
## to character vector
## Warning: Column `names`/`commcode` joining factor and character vector,
## coercing into character vector

## Warning: Column `names`/`commcode` joining factor and character vector,
## coercing into character vector

## Warning: Column `names`/`commcode` joining factor and character vector,
## coercing into character vector

Modeling options

Districts with & without CTAR

A Random Effect by Catchment

Models with a random intercept by catchment. May help to capture variation between clinics not explained by travel times alone. And provide more realistic estimates of vials.

Model Results

Convergence fig

DIC fig

Estimates

Predictions to Data

Predictions out-of-fit

Moramanga estimates on Mada commune and districts

Mada commune and districts estimates on Moramanga data

Burden estimates